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A GROWTH MODEL IN A RANDOM ENVIRONMENT 



Janko Gravner, Craig A. Tracy, Harold Widom 



1. Introduction. 

Processes of random growth and deposition have a long history in the physics hterature, 
typically as models of systems far from equilibrium (e.g., [Mea] and the more than 1300 refer- 
ences listed therein) . They made their appearance in probabilistic research about 35 years ago, 
with arguably the most basic growth rule, first passage percolation ([HW]). The fundamental 
asymptotic result is an ergodic theorem: scaled by time t, the growing set of sites approaches a 
deterministic limiting shape. As these early successes were based on nonconstructivc subaddi- 
tivity arguments, they posed two natural questions: (1) can the asymptotic shape be identified 
analytically and (2) how large are fluctuations about the limit? While there has been no res- 
olution of the first issue, ingenious probabilistic and geometric arguments have yielded much 
progress on the second ([Ale]), although the matter is still far from settled. It is therefore of 
some importance to be able to provide a complete answer on some other simple, but nontrivial, 
interacting growth process. It turns out that several two-dimensional oriented models with last 
passage property ([Sepl], [Sep2], [Johl], [Joh2], [BR], [PSl], [GTWl]) are most convenient, as 
they can be represented, on the one hand, as particle systems related to asymmetric exclusion 
and, on the other hand, as increasing paths in random matrices and associated Young diagrams. 
This allows explicit answers to both questions (1) and (2) above. 

In this paper we continue to study Oriented Digital Boiling (ODB) (Feb. 12, 1996, Recipe at 
[Gri], [Gra], [GTWl]), perhaps one of the simplest models for a coherent growing interface in 
the two-dimensional lattice 7?. The occupied set, which changes in discrete time t = 0, 1, 2, . . . , 
is given hy At = {{x,y) : x G Z, y < /it(x)}, and the height function ht evolves according to the 
following rule: 



Here e^^t are independent Bernoulli random variables, with P{£x,t = 1) = Px- Thus the prob- 
ability of a random increase depends on the spatial location. It remains to specify the initial 
state which will be 



In [GTWl] we analyzed the homogeneous case Px = p, identifying the following four asymp- 
totic regimes: 



ht+i{x) = max{/it(x- l),ht{x) + £x,t}- 



(1.1) 
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- Finite x GUE Regime: if x is fixed and t — oo, then {ht{x) — pt)/y/p{l — p)t — >■ M^, 
a Brownian functional whose law can be computed explicitly as the largest eigenvalue of 
(x + 1) X (x + 1) hermitian matrix from the Gaussian Unitary Ensemble (GUE). 

- GUE Universal Regime: if a; is a positive multiple of t, and x/t < 1 — p, then there exist 
constants ci and C2 so that {ht{x) — Cit)/{c2t^^^) converges weakly to a distribution F2 
known from random matrix theory ([TWl]). 

- Critical Regime: if x = (1 — p)t + o{^/i), then P{ht{x) — {t — x) < —k) converges to a 
k X k determinant. 

- Deterministic Regime: if x is a positive multiple of t, and x/t > l—p, P(ht(x) = t — x) — > 
1 exponentially fast. 

The focus of this paper is ODB in a random environment, in which p^ are initially chosen at 
random, with common distribution given by P{px < s) = F{s). We will also assume that p^ are 
independent, although in several instances this assumption can be considerably weakened. In 
statistical physics, processes in a random environment are often called disordered systems, or, es- 
pecially in the Ising-type models, spin glasses. In this context, the random environment (choice 
of Px) is referred to as quenched randomness, as opposed to the dynamic {thermal) fluctuations 
induced by the coin flips ex,t- In general, rigorous research in this area has been a notoriously 
difficult enterprise; for some recent breakthroughs (as well as reviews of the literature) we refer 
the reader to [SK], [NS], [NV] and [Tal]. 

We now state our main results. Throughout, we will denote by ( • ) integration with respect 
to dF and p a generic random variable with distribution F. 

Construct a random m x n matrix A = A{F), with independent Bernoulli entries Sij and 
such that P{sij = 1) = pj, where, again, pj = p are i.i.d. Label columns as usual, but rows 
started at the bottom. We call a sequence of I's in A whose positions have column index 
nondecreasing and row index strictly increasing an increasing pa,th in A. Let H = H{m,,n) be 
the length of the longest increasing path. (Sometimes, to emphasize dependence on F, we write 
H = H[F) = H[m, n, F).) The following lemma is then easy to prove ([GTWl]). 

Lemmal.l. Under a simple coupling, ht{x) = H{t — x,x + 1). 

We will therefore concentrate our attention on the random matrix A from now on, switching 
to the height function only occasionally to interpret the results. We also note that Lemma 1.1 
demonstrates that ODB is equivalent to the Seppalainen- Johansson model ([Sep2], [Joh2]). 
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Our first theorem identifies the time constant. In the sequel, we will present two completely 
different methods for proving these limits, a variational approach and a determinantal approach. 
The first method (which is similar to the one in [DZ]) is based on the crucial symmetry property 
of H (Lemma 2.2) and provides some information on the longest increasing path itself, while the 
second one is deeper and more precise and thus able also to determine fluctuations. The paper 
[SK] studies a related model, presents yet another technique, based on an exclusion process 
representation, and observes similar phase transitions. Throughout this paper, we let 



be the right edge of the support of dF and assume that n = am for some < a < cxd. (Actually, 
n = [am\ , but we drop the integer part as it obvious where it should be used and to avoid 
complicating expressions.) We also define the following critical values 



Note that c determines the limiting shape of At, namely lira At /t, as t ^ cx) for the corner 
initialization given by (1.1). By virtue of the Wulff transform, it then also gives the speeds of 
some half-planes, i.e., lim At /t when comprises points below a fixed line. See [SK] for much 
more on this issue. 

Theorem 1. The limit in (1-3) exists almost surely. Ifb=l, then c{a,F) = 1 for all a, while 
if b < 1, then 



b = b{F) = mm{s : F{s) = 1} 



(1.2) 




and define c = c{a,F) to be the time constant 



(1.3) 



c = c{a,F)= lim — . 

m— >oo m 





Note that that (^{b — p) ^) = oo iff = iff there is only one critical value. 
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Next we turn our attention to fluctuations. In this paper we present complete results for the 
pure regime a!^ < a < and for the (easy) deterministic regime < a. The composite regime 
a < is addressed in [GTW2], while both critical cases when a equals either critical value 
currently remain unresolved. To explain the results, and connect with the spinglass terminology 
we have just used, we turn to a simulation. For an example, we use F{s) = 1 — (1 — 2s)'^ so 
that b = 1/2, ac ~ 6.3 and a'^ ~ 0.5 and run the simulation until time t = 40, 000 (with a single 
realization of the environment and the coin flips). When x is close to the origin, it is clear from 
the picture that the interface mostly consists of sheer walls followed by flat pieces. The walls 
correspond to the rare sites with update probability Px close to 1/2. Those are much faster than 
the other sites so they pull ahead of their left neighbors, creating walls, and dominate their right 
neighbors by "feeding" them at nearly largest possible rate. In fact, this state of affairs persists 
up to about X = t/3 although close to x = t/3 these effects are less pronounced. In the pure 
regime, when x/t ranges approximately from 0.333 to approximately 0.863, the fluctuations are 
much more regular, and in fact, as we will demonstrate, asymptotically normal. For larger x/t 
the shape has slope —1 and no fluctuations. 



25000 




Figure 1. Two ODB simulations, as explained in the text. 



For comparison, consider the case when p is uniform on [0,1/2], the case that has Oc = 
l/(ln4 — 1) ~ 2.59 and a'^ = 0. The fluctuations arc normal up to z/t ~ 0.72. Figure 1 depicts 
the results of simulations, first complete boundaries of two occupied sets (the top curve is the 
uniform case), then two details (the right curve is the uniform case) for x G [1000, 5000]. 



5 



Theorem 2. Assume that b < 1 and oi^. < a < a^- Let a be as in Theorem 1 and 



Then, as m ^ oo, 



H — cm 
T^/a ■ m^/^ 



iV(0,l). 



Assume that p is uniform [0, 1/2] to illustrate Theorems 1 and 2. Together they imply that 
there exist ci and C2 so that {ht{x) — cit)/{c2t)^^'^ — ^ -^(0, 1), where ci determines the limiting 
shape and C2 is the variance. These two quantities are presented in Figure 2, ci is the top and 
C2 is the bottom curve. For comparison, the shape of homogeneous ODB with = {p) = 1/4 is 
also drawn (middle curve). Note that ci and C2 approach 1/2 and 1/4, respectively, as a — 0, 
indicating that for small x/t the interface growth is governed by the largest update probability, 
which is close to 1/2. Finally, we do the same computation for the other example in Figure 1. 
The variance is now drawn only on [a^, 1]. 




Figure 2. ci (top), C2 (bottom) and the shape for p^ = {p) (middle) vs. x/t. 
The two distributions are uniform [0,1/2] (left) and F{s) = 1 - (1 - 2sf . 



We note that both a.s. convergence to the limiting shape (which is equivalent to a.s. con- 
vergence in (1.3)) and its convexity follow from subadditivity, which in turn is a consequence 
of the fact that this is an oriented model in which influences only travel in one direction. To 
be more precise, fix integer sites (xi,?/i), (a;2,y2) G Z_|_ x Z_|_ and define times T'(xi,yi),(x2,y2) 
follows. First wait until time T'(o,o),(xi,yi) when the dynamics reaches (xi,yi). Then restart the 
dynamics from the initial state 



hci{x) 



yi, ifa; = xi, 
— oo, otherwise. 



and let T/ 



(xi,yi),{x2,y2 



) be the time at which the occupied set reaches {x2,y2)- This random 



< T, 



(0,0),(a;i,{/i) 



+ T, 



(xi,yi),{x2,y2)- 



variable is independent of Px for x < xi — 1 and ^(o,o), (0:2,1/2) 
Therefore, the subadditive ergodic theorem can be applied as in the first chapter of [Dur]. 
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The main step in the proof of Theorem 2 estabhshes a hmit law for fluctuations conditioned 
on the state of the environment. In many ways, such a result is more pertinent to understanding 
physical processes modeled by simple growth models such as ODB. 

Theorem 3. Assume that b < 1 and o!^ < a. < a^- Then there exists a sequence of random 
variables Gn G cr{pi, . . . and a constant 5o 7^ (both depending on a) such that, asm ^ oo, 

P{ ^-1 'iZ < s \ Pi,...,Pn] F2{s), 
almost surely, for any fixed s. 



The random variables Gn = Cnm are given in terms of the solution of an algebraic equation 
in which pi, . . . ,Pn appear as parameters (see (3.4) and (3.5)), while the deterministic constant 
go is specified before the statement of Lemma 3.5. The limiting distribution function F2 first 
arose in connection with eigenvalues of random matrices ([TWl], see [TW2] for a review). Since 
then it has been observed in many other contexts, including growth processes ([Johl], [Joh2], 
[BR], [GTWl], [PSl], [PS2]). Most suitable for computations is the identity 

F2{s) = exp J ~ s)q{x)^ dx 
where q is the unique solution of the Painleve II equation 

q" = sq + 2q\ 

which is asymptotic to the Airy function, q{s) ~ Ai(s) as s ^ 00. When proving limit laws, 
it is more useful that F2 can be represented as a Fredholm determinant (see e.g. [GTWl] and 
Section 3 below). 

In Theorem 3 the environment is assumed as given, H is approximated by the quenched shape 
Gn, with the fluctuations about this shape of the order m^/^ and given by the F2 distribution. 

As we prove in Section 3, {am)-^/'^{Gr, — cm) converges to the standard normal, making it 
clear why Theorem 2 holds: the environmental noise eventually drowns out the more interesting 
quenched fluctuations of Theorem 3. An illustration is provided in Figure 3, in which p is 
again uniform on [0, 1/2] and ht (solid curve), deterministic approximation based on Theorem 2 
(dotted curve), and the much better random approximation based on Theorem 3 (dashed curve) 
are all depicted at times t = 100, 200, . . . , 1000. 
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Figure 3. Approximations to ht (solid curve) based on 
Theorems 2 (dotted curve) and 3 (dashed curve). 



In conclusion, we note that the connection between random matrix theory and random com- 
binatorial objects, which has become the key to rigorous understanding of random interface 
fluctuations, made its initial appearance in [BDJ], while an inhomogeneous model of ODB type 
was first studied in [ITWl]. This last paper, together with its companion [ITW2], extends the 
study of random words from the homogeneous case in [TW3] in a somewhat analogous way 
as the present paper builds on the work in [GTWl]. In particular, connections with operator 
determinants (from the beginning of Section 3) arc very similar (see also [Rai] which features a 
general inhomogeneous setup). However, randomness of the environment, which seems to be a 
new feature in rigorous analysis of explicitly solvable models, then forces our techniques to take 
a novel turn. 

2. A variational characterization of the time constant. 

We start by a remark on constructing the random matrix A. The most convenient design uses 
as the probability space (O, P) a countably infinite product of unit intervals [0, 1] with Lebesgue 
measure. A copy of the unit interval (and thus a factor in the product) is associated with each 
point in N X N and in addition, with each positive integer in N. (The former factors correspond 
to matrix entries, and the latter to its columns.) If a; = {mij,Cj) G ^2 is a generic realization, 
we define the following random variables: pj = F~^{cj) (where F~^(x) = snp{y : F{y) < x} as 
usual) and Sij = l^mijKpj}- By restricting to the m x n rectangle at the lower right corner of 
N X N, this constructs the random matrices A for all m and n simultaneously. The following 
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useful lemma also follows immediately. 

Lemma 2.1. If Fi < F2 are two distribution functions, the two corresponding random matrices 
A{Fi) and A{F2) can be coupled so that H{F2) < H{Fi). 

Next we state the crucial property for the variational approach to work: conditioned on the 
environment, H is a symmetric function of flip probabilities. 

Lemma 2.2. A regular conditional distribution 

P{H <h\pi,...,p^) 

is a symmetric function of pi, . . . Pn- 

Proof See section 2.2 of [GTWl]. □ 

Somewhat loosely, wc denote by Hn the random variable H obtained by fixing pi, . . . ,pn- In 
fact this is nothing more that a shorthand notation, e.g., E{ip{Hn)) = E{ip{H) \ pi, . . . ,pn) for 
any bounded measurable function tp. 

The time constant c{a, x) = c(a, S^) for the case pj = x is given in [GTWl]. The next lemma 
summarizes the relevant conclusions. 



Lemma 2.3. Assume that dF = S^- Then 



c = c{a, x) 



1\foL\Jx{\ — a;) + (1 — a)a;, (1 — a;)/a;>a, 
1, (1 — x)/x < a. 

Moreover, for every £ > there exists a constant 7 = 7(e) > so that 



(2.1) P{\H/m-c\> e) Ke-^"" 

for m > mo(£, a, x). 

Proof. The formula for c follows from (3.1) in [GTWl], while the large deviation estimate can 
be proved by the method of bounded differences as in Lemma 5.4 of [Gra]. □ 



It turns out the following function is more convenient than c. 
C,{y,x) = yc{l/y,x) 



( 2^a/x(1 - x) + {y- l)x, x/ (l-x) < y, 
I y, x/{l - x)>y. 
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Note that the partial derivative 

( y-'^/'^^/xiT^^ + x, x/{l-x)<y, 
I 1, x/{l -x)>y. 

is decreasing in y (obviously) and increasing in x (easily checked). In particular, C(")^) is a 
convex function. 

We now derive a variational problem for c, initially without paying attention to rigor. Start 
by a nice distribution function F and approximate it by the discrete distribution function given 

by 

F(p,4)=Aft(.) = F(')-F(i^),i = l....... 

Let ip : [0, a] — [0,1], ?/'(0) = 0, ip{oi) = 1 be a nondecreasing function, with A'0fe(z) = 
ip{aF{i/k)) — ip{aF{(i — l)/k)). Define the functionals: 



jr(V;)= / C{'ip'{aF{x)),x) -adFia 
Jo 



and 

k 



\aAFk{i)'k 

1=1 ^ ^ ' 

Generate the pj's and denote by Ni the number of pj equal to i/k. By Lemma 2.2, we 
can assume flip probability 1/A; in the first Ni columns, 2/k in the next N2 columns, etc. 
Moreover, the strong law suggests that the identity Ni = AFk{i)n nearly holds. As we know the 
asymptotics for the longest increasing paths in the slivers of widths iVj in which the probabilities 
are constant, the longest increasing path in A is determined by the most advantageous choice 
of transition points between the slivers. These transition points arc specified by a function tp as 
described above. If we approximate the differences with derivatives, we obtain 

c(a,F) = lim c{a,Fk) 

fc— >oo 
k 



lim max > c — — — -— , — • /\wu(i) 



= lim meixJ-k{ip) 



= Jim ma^^C ('ip'{oiF{i/k)),^^ -aF'ii/k)-^ 

^ i=l ^ ^ 
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At this point, we remark that a connection between longest increasing paths and variational 
problems has appeared before in the literature. The result closest to ours is by Deuschel and 
Zeitouni ([DZ]), who used a variational approach to study a variant of Ulam's problem. In 
their case, a number of points in the unit square is chosen independently according to some 
distribution with a density, then a longest sequence, increasing in both coordinates, is extracted 
from this sample. The Deutehel-Zeitouni functional is different from ours as the length of the 
longest increasing path has a non-trivial dependence on a (that is, through c) in our case. 

The (integrated) Euler functional for the variational problem is 

Cyii^'{x),F-\x/a)) = a, 

or, writing g{x) = ip' {aF{x)), 

(2.2) Cy{g{x),x)=a. 

Since Cy < 1 and equal to 1 if and only if x/(l — z) > y, the integration constant a G [0, 1]. If 
o = 1, then g(x) <x/(l — x), ({g{x),x) = g{x) and 



c{a,F)=[ 'ilj'{aF{x))adF{x) = 1. 
Jo 



Assume now that 6 < 1. In this case, it is necessary to specify g only on [0, b). But (2.2) gives 

The constant a is given by the boundary conditions. Assuming that (2.3) holds on [0,6], 

(2.4) l = a [ g{x)dF{x) =a [ ~ dF{x). 

Jo Jo {(^ ~ ^) 



The smallest the last integral can be is when a = 1, which yields the condition 

rb 

X 

^ X ' ' ar 



l>a [ dF(x) = " 

Jo 



On the other hand, the largest that the integral in (2.4) can be is when a = b. Therefore, if 
a G (a^, Oc)) we have found the minimizer and 



c{a,F)= / (^{g{x),x)adF{x) = a 
Jo 



— a^p + 2ap 
(a - p)2 



which reduces, upon using the defining equation for a, to the formula in Theorem 1. 
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lia < a'f., the minimizer ip has to make a jump of size 1—a/a'^at a. The natural interpretation 
for this is that the minimizer given by (2.3) is used in the lower left part of A with dimensions 
{a/a'^)m x (n — 1). To the resulting increasing path in this submatrix one needs to add the 
number of I's in the upper segment of length (1 — a/a'^)m in the last column, in which nearly 
the largest probability b is used. Therefore, 

c{a,F) = c{a',,F).^+b(l-^), 

which again reduces to the appropriate formula in Theorem 1. 

We now proceed to give a proof Theorem 1, the heart of which is a somewhat involved 
multistage approximation scheme. 



Proof of Theorem 1 when 6=1. This follows simply by observing that, for any e > 0, maxj pj 
b a.s. as m — oo. Since a trivial lower bound is obtained by using only the column with the 
largest pj, one concludes that liminf iJ/m > b a.s. □ 



Proof of Theorem 1 when a G (a^, ccc)- 
We begin by the following lemma. 



Lemma 2.4. Assume that a sequence of distribution functions Fn converges to F in the usual 
sense (i.e., the induced mea,sures converge weakly). Assume also that 6(-Fjv) — b{-^) ^'"'^ ^^^^ 
aN — «• Then c(aAr,-Fjv) c(a,F) (as given in Theorem 1). 

Proof. If a' > b{F) and 

„/x(l-x)(.'-x)-.FW>l. 

then for a large A^, a' > b{FN) and, since the integrand is bounded, 

aN J x{l — x)(a' — x)~^ dFN{x) > 1. 

Hence On = a{aN, F^) > a' . If ooi then x{l—x){aN—x)~'^ converges to x(l— a;)(ao— a;)~^ 

uniformly for x € [0, a'] and so 

1 = ajv J x{l — x)(ajv — x)~^ dF^^x) a J x{l — x){ao — x)~^ dF{x). 
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Therefore oq = a{a,F) and consequently qn — a{a,F). As x{aN — x) ^ also converges uni- 
formly on [0, a'], 

c(q;jv, -Pjv) = Oat + CKjv(l — Oat) J x{aN — x)~^ dF^ 

^a + „(l-»)/x(a-.)-'<iF = c(„,n 

□ 

First we assume that F is nice, that is, a one-to-one function on [P, b] C (0, 1), with F{P) = 0, 
F{b) = 1, and continuously differentiable on (0, 1). We also assume that ^ is the class of non- 
decreasing convex functions ip G C^[0, a], with ^(0) = 0, 4>{a) = 1, tp'i^) ^ P/^- This last 
assumption is necessary because C,{y-,x) is not Lipshitz near y = 0. 

Lemma 2.5. Assume that a G {a,'^,ac)- Among all ^ G the functional is uniquely 

maximimized by 

^{x)= f g{F-\u/a)fdu, 
Jo 

where g is given by (2.3). 

Proof. This follows from standard calculus of variations. Both 'tp'{0) > /3/2 and convexity of ip 
are easily checked. □ 

We now justify the approximation steps in the heuristic argument, using the same notation. 
First, if e > is fixed, then with probability exponentially (in n) close to 1, 

(1 - e)AFk{i)n < Ni < {1 + e)AFk{i)n 

for every i = 1, . . . ,k. By obvious monotonicity, the longest increasing path in A is then bounded 
above by the longest increasing path in A' in which all iVj = (1 -|- £)AFf.{i)n, and therefore we 
can get an upper bound by increasing a to a(l -|- 2e) and assuming iVj = AFk{i)n. A lower 
bound is obtained similarly. As our final characterization of c is continuous with respect to a 
(Lemma 2.4), we can, and will, assume that iVj = AFf.{i)n from now on. 

The above paragraph eliminates randomness of Pj's; we now proceed to replace the coin 
flips with deterministic quantities. Again, fix an e > and let M = em. For ji < j2 and 
i = 1, . . . , n, consider the longest increasing paths vr^j ja^j between {Fk{i — l)n, ji) (noninclusive) 



13 



and {Fk{i)n,j2) (inclusive). Then, with probabihty exponentially close to 1, the length of any 



^ji,i2,i is most 



n. 



(This uses Lemma 2.3 when j2 — ji is divisible by M and fills the rest by monotonicity. Note 
that Lemma 2.3 is therefore only applied finitely many times for fixed e and A;.) The lower 
bound is obtained by rounding down instead of up. It follows that the length of any T^j^,j^,i is 
bounded above (resp. below) by 

(2-5) cffexrVAF^W"- 



AFfc(i)n 

computed on the matrix of size (m + M) x n (rcsp. (m — M) x n). Once again we can use 
continuity to assume that the length of any nj^^j^^i is given by (2.5). 

It remains to show that the discrete deterministic optimization problem max^ (V") is for 
large k close to its continuous counterpart max^ ■?^(V')- To this end, we first prove that we can 
indeed restrict the set of function ip to those in ^, i.e., those that are convex and have a large 
enough derivative. Let Axi = aAFfc(i), Ax2 = aAFf~{i + l), Ayi = Atl;f~(i), Ay^ = A'i/;fc(i + l), 
Ay = Ayi + Ay2, pi = i/n, P2 = {i + l)/w. Then 

M ,(||i.,,)^.„c(^..)A. 

is nondecreasing with decreasing Ayi as soon as pi < p2 and Ayi/Axi > Ay2/Ax2. This means 
that the maximum is achieved at a convex tp. Similarly, the expression (2.6) is nondecreasing 
with increasing Ayi if pi > (3 and Ayi/Aa;i < 5/(1 — 6), and therefore the maximum is achieved 
at aip e^. 

Next we note that 

k 

^fe(V') <Y,ai^'iaFii/k)),i/k)aAFk{i), 

i=l 

while 

k 

^(0) >J2Ci^'i^Hi/k)),i/k)aAFk+,(i). 

i=l 

Therefore, max^ Tki'tp) < max^ J^{ip) + 0(1/ k). As a lower bound is obtained similarly, this 
concludes the proof for nice distribution functions F. 

To prove the general case, we again use Lemmas 2.1 and 2.4. For an arbitrary distribution 
function, choose nice F^ so that F^ < F and F-l(i/jv,i] 

< F+ and F^ ^F and 6(F^) ^ 6(F). 
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Then c(q;, Fjy) — c{a,F). By Lemma 2.1, it immediately follows that lim sup i?/m < c{a,F) 
a.s. 

The lower bound, however, does not immediately follow as F is not below F^j. The remedy for 
this is to assume that F{l/N) < 1/2, replace a with a' < a, and observe that the distribution F 
will induce, with probability exponentially close to 1, at least {a — a')m/4 probabilities pj > 1/A^. 
Therefore the length of the longest increasing path in a m x am matrix using F is eventually 
above the length of the longest increasing path in a m x a'm matrix using F^. By Lemma 2.4, 
liminf ii"/m > c{a,F) a.s. □ 



Proof of Theorem 1 when a < a'^. 

Applying the same strategy as before we construct sequences {F^} of distribution functions 
which satisfy F^ < F < F^ and for which Theorem 1 already holds, and such that c(F^) 
and c(-P^) approach the same limit as N oo. Lemma 2.1 will then complete the proof. (We 
suppress a from the notation, since it is the same throughout this proof.) 

Take a sequence r?7v \ such that b — r]N are points of continuity of F. Let F^ agree with 
F outside [b — r)N,b), while on [b — r]N,b) the two functions are constant: F^ = F{b — rjN) and 
F^ = 1. Let Eat = 1 — F{b — tjn)', note that Eat — and dF^ = l(o,6-r;jv) + ^N^b and 
dF^ = 1(0, b-riN) + Sb-riN- Clearly the already proved part of Theorem 1 applies to both 
F^ and F^. 

We proceed to show that a{F^) b. If this does not hold, the fact that a{F^) > b{F^) = b 
implies that there exists an 77 > so that a{F^) > b + tj along a subsequence. Then S = 
{p{l -p) [ (6 - p)'"^ - {b + rj- p)"2 ] ) > and 



(2.7) l = a [ 

Jo 



^-"^ .x(l-x) ,^ 6(1-6) 
^ ^ dF + asN 



b 



{a{F^)-xY {a{F^)-h) 
x{l-x) , _ 6(1-6) 



,2 



^ « / IT-, ^ dF + a£N 

Jo {b + V- xY V 

, a 6(1-6) 
< -0 + — + aeN 5-^, 

along the same subsequence. As N ^ 00, this yields a contradiction with a < a!^. 
Now 

(2.8) c{F^) = a{F^) + ail-a{F^)){ / — — dF + SN 



a{F, 



N > 
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By (2.7), 



To show that 

(2.9) (l|p<,_,„}p/(a(F^) -p)) ^ {p/{b-p)) 

we note that the integrand on the left of (2.9) is uniformly integrable (as it is bounded by 
p/ ~ P)i which is square-integrable) and converges to the integrand on the right a.s. By (2.8) 
and (2.9), 

c{F^)^h + a{l-h){p/{l-p)). 
The argument for c{F^) is very similar and hence omitted. □ 



Proof of Theorem 1 when a > ac- If a /■ ac, then a{a, F) y 1 and hence c{a, F) y 1. D 



We note that the above proof of Theorem 1 actually shows exponential convergence to c, 
that is, (2.1) in Lemma 2.3 holds in random environment as well. Also, once probabilities are 
ordered using Lemma 2.1, one could investigate convergence, in the sense of [DZ] and [Scpl], of 
a longest increasing path in A to the maximizer of ^(V') . This is easy to prove if is nice (cf. 
Lemma 2.5), but it actually holds whenever the maximizer is unique. 

We conclude this section by showing that the deterministic case indeed has no fluctuations. 



Proposition 2.6. Assume that b < 1 and a > ac- Then P{H = m) converges to 1 exponentially 
fast (and therefore P{H = m eventually) = IJ. 



Proof. We begin by modifying the construction from Section 3.3.1 of [GTWl]. Recall that 
random mxn matrix is the lower left corner of an infinite random matrix. For an G N x N, 

let ?7(i,j) = infjfc > 1 : £(j_|_fc^j) = 0} be the relative position of the first above and 
~ ^ 1 • ^(i,j+k) = 1} the relative position of the first 1 to the right of 

Now define i.i.d. two-dimensional random vectors Xi = (^i, 771), X2 = (42,772), ... as follows: 
6 = 4(0,1)' Vi = '7(6,1)' 

6 = 4(^i,i+,7i), m = V{^i+i2,i+m)^ 

^3 = 4($i+«2,i+»7i+';2)' ^2 = r]{ii+i2+i3A+vi+m)' 
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Let Sk = (0, 1)+Xi H \-Xk be the corresponding random walk, and T^, (resp. T^) be the first 

time Sk is in {{x,y) : x > n} (resp. {{x,y) : y > m}). If < r„ then there is an increasing 
path of I's inside the m x n rectangle which goes through its "roof" without skipping a row, 
thus 

{H<m}c {r„ < T;,}. 

Therefore, we need to show that P(T„ < T^) goes to exponentially fast. To this end, note 
that, for any £ > 0, 

oo 

(2.10) P(r„ < T;,) < P{Tn A r4 < em) + ^ P(6 + • • • + ^ > a{Vi + ■■■ + %))• 

k=sm 

If we show that and 771 have exponential tails, and that -B(^i) — aE{r]i) < 0, the we can 
choose a small enough e > so that the upper bound in (2.10) decays exponentially. First, 
-P(Ci > k) = {1 — p)^~^ and so -E(^i) = 1/ {p)- Moreover, the conditional distribution of p given 
that a single coin flip gives 1 is 

dFi{x) = -j-^xdF{x), 

therefore 



P{m >k)=[ x''-^dFi{x) 
Jo 



and so E{r]i) = {p/{l - p)) / {p). □ 



3. The saddle point method and fluctuations. 

Throughout this section, we assume that 6 < 1 and that a = n/m is fixed (but see Remark 
3 at the end). In addition, our standing assumption will be that 

a'^ < a < OLc- 

We will investigate the limiting behavior of P{H. < h) without using results proved in Section 
2. An asymptotic analysis of this quantity when a < is carried out in [GTW2]. 

Wc begin with deterministic inhomogcncous ODB, in which the jth column is assigned a 
fixed deterministic probability pj. At first, our derivation will use a fixed n and no particular 
properties of the eventual random choice of the environment. For notational convenience, we 
therefore drop the subscript n, which practically every quantity would otherwise have. See the 
discussion preceding the key formula (3.6), where the random environment is reintroduced. 

As explained in [GTWl], Sec. 2.2, we have 

p{H<h) = \{{i-PirDM, 
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where is the hx h ToepHtz determinant with symbol 

n 
J = l 

and rj = pj/{l—pj). Applying an identity of Borodin and Okounkov ([BO], see also [BW]) this 
becomes 

P{H <h) = det (/ - Kh), 
where is the infinite matrix acting on ^^(Z"*") with j, k entry 

oo 

Kh{j,k) = ^(</j_/<<p+)/,+j+^+i {^^/^-)-h-k-e-i- 

The subscripts here denote Fourier coefficients and the functions (p± are the Wiener-Hopf factors 
of if, so 

n 

^+iz) = l[il + rjz), (^_(z) = (l-z-^) — 
The matrix is the product of two matrices, with j, k entries give by 

(y,+/<^_)_,_,_fc_i = ±-Jl[{l + r,z) {z - ir dz 



and 



= ^/ \{{l + r,z)-^ {z - \)-^ z^-^-i-^--" dz. 



The contours for both integrals go around the origin once counterclockwise; in the second integral 
1 is on the inside and all the —t^^ are on the outside. 

Eventually we let m, n ^ oo and will take h = cm + sm}-!^ where c, as yet to be determined, 
gives the transition between the limiting probability being zero and the limiting probability 
being one. In [GTWl] we considered the case where all the pj were the same. We found that 
with c chosen as in Lemma 2.3 we could do a steepest descent analysis. The conclusion was 
that the product of the two matrices scaled, by means of the scaling j —>■ m^^^x, k —>■ m^^^y, 
to the square of the integral operator on (0, oo) with kernel Ai{gs + x + y), where g is another 
explicitly determined constant. This gave the limiting result 

lim P{H <cm + sm^^^) = F2{gs), 

n—>oo 

where i*2(s) is the Predholm determinant of the Airy kernel on (s, oo). We can do very much 
the same here, li h = cm + sm}/^ and we set 



'^{z) = \{{l + r,z) [z-l] 
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then 



(3.1) {ip^/^.).h-j-k-i = ^.j ^{z)z^^'''+^+'^dz, 



(3-2) (<^_/(/j+)^+j+fe+i = y V'W ^ 

To do an eventual steepest descent we define 

1 a " 

a{z) = — log il){z) = — ^ log (1 + rj,2) + log - 1) + (c - 1) log z, 
m n .^^ 

and look for zeros of 

3.3 a'z=-7 — + 7 + ■ 

^ ' ^ ' n l + r^z z-\ z 

The number of zeros equals one plus the number of distinct rj. There is a zero between two 
consecutive l/rj and, in general, two other zeros which which are either unequal reals or a pair 
of complex conjugates. In the exceptional case there is a single real zero of multiplicity two. We 
choose c so that we are in this exceptional case. If the double zero is at z = ti then u and c must 
satisfy the pair of equations 



a V— V r, 1 c — 1 „ 
— > 1 1 = 0, 

j=i ■> 

n / X 2 -I 

\ + ^ + 



^ = 0. 



3- 

If we multiply the second equation by u and subtract we get 



(3-4) - E TT^ ' 

^ ^ n ^ (l + r-i 



(1 + rjuY {u - If 
The first equation gives 

(3.5) c = ^-^y- 

^ ' l-u n^l 



Conversely, if the second pair of equations is satisfied then so is the first. 
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Lemma 3.1. Assume that an ^ < 1 and set u = max{— l/rj}. Then equation (3.4) has 
a unique solution u G (u, 0) and if c is then defined by (3.5) we have c G (0, 1). 

Proof. The left side of (3.4) decreases from oo to an~^ Yl^j ^ ^ runs over the interval {u, 0) 
whereas the right side increases and has the value 1 at = 0. Our assumption implies the first 
statement of the lemma. As for the second, c > since u < and each 1 + rjU > 0. Moreover, 
Schwarz's inequality, our assumption and (3.4) give 

l + rjU~ ^ ) {1 + VjuY J l-u 

Hence 

1 u 

c<- = 1. 

1 — u 1 — u 

□ 

To derive the asymptotics using steepest descent we have to compute a"'{u) and understand 
the steepest descent curves. For the first we multiply (3.3) by differentiate twice and use the 
fact that cf'{u) = (t"{u) = to obtain 

"'( ^ 2a " r| 2 

Note that a'" (a) > since u<0. 

There are three curves emanating from z = u on each of which SJcj is constant. One is 
= 0, which is of no interest. The other two come into u at angles ±7r/3 and ib27r/3. Call the 
former C"*" and the latter C~ . Approximate shapes of these curves are illustrated in Figure 4. 
For the integral involving ip{z) we want |V'(-2)| to have a maximum at the point u on the curve 
and for the integral involving 'ip{z)~^ we want |V'(-2)| to have a minimum at u. Since a"'{u) > 
the curve for ip{z) must be C"*" and the curve for ^[z)~^ must be C~ . 

As for the global natures of the curves, can only end at a zero of 'tp{z)'^^, at a zero of 
o-'{z), or at infinity. The two curves are simple and cannot intersect since \ip{z)\ is decreasing 
on as we move away from z = u while \ip{z)\ is increasing on C~. It follows that C+ closes 
at z = 1, while the two branches of C~ go to infinity. Prom the fact that 

n 

9 a{z) = — arg (1 + rjz) + arg (z — 1) + {c— 1) arg z 
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is constant on C~ we can see that the two branches go to infinity in the directions argz = 
■^ci^l{c + a{\ — u)), where u is the fraction of rj equal to zero (which is the same as fraction of 
the pj equal to zero). Observe that in the integral in (3.1) the path can be deformed into (7+ 
and in the integral in (3.2) the path can be deformed into C~ . Both contours will be described 
downward near u. 




u 1 



Figure 4. The steepest descent curves as described in the text. 



To see formally what steepest descent gives, we replace our matrices M(j, k) depending on 
the parameter m and acting on £^(2+) by kernels m^/^M{m^^^x, m}^^y) acting on L^(0, oo). 
Thus (3.1) becomes the operator with kernel 



2m J 



If steepest descent worked, the main contribution would come from the immediate neighbor- 
hood of z = u. We would set z = u + (,, make the replacements 

aiz)^a{u) + ^a"\u)e, z^ue^'^, 
in the integral and integrate (downwards) on the rays arg C = ib7r/3. The above integral becomes 

and we can then replace the rays by the imaginary axis (downwards). The variable change 
C — —iC,/m}/^ replaces this by^ 

_^ma{u) ^m^^^x / ^^cr'" (u)C^ -i{s+x+y)C/u 



iRecall that the Airy function is defined by Ai(a;) = ^ 
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= _gma(«) ^m'/^(s+x+y) \u\g Ai{g{s + X + y)), 

where we have set 



g = \ l^"'{u) 



-1/3 



Thus, if we multiply the matrix entries on the left side of (3.1) by 



then the result has the scaling limit the operator on L^(0, oo) with kernel 

\u\gA\{g{s + x + y)). 
Similarly if we multiply the matrix entries on the left side of (3.2) by 

then the result has the scaling limit the operator on -L^(0, oo) with kernel 

\u\~^gA\{g{s + x + y)). 

It follows that the product of the two matrices has in the limit the same Fredholm determinant 
as the operator with kernel 

\ K\{g{s + x + z)) Ai(5(s + z + y)) dz 
Jo 

nOO 

= g Ai{g{s + x) + z) ki{g{s + y)) dz 
Jo 

which in turn has the same Fredholm determinant as the kernel 

roo 

/ Ai{gs + x + z) Ai(gs + z + y) dz. 
Jo 

This Fredholm determinant equals F2{gs). 

Assuming the argument we sketched above goes through we will have shown that, in some 
sense, 

lim P{H < cm + srn^^^) = F2{gs), 

n — >cc 

where c and g are as above and determined once we know the pj and a. 

We begin the rigorous justification by introducing some notation. Recall that we consider a 
random environment in which the probabilities pj are chosen independently with distribution 
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function F. We explained the notation iJ„ after Lemma 2.2; in addition, we give the subscript n 
to the quantities (T„(2;), u^, gn-> curves to emphasize that they are functions of pi, . . . 
Therefore 

P{H <h) = {P{H^ < h)) , 
where ( • ) is the expected value with respect to pi, . . . 

Our object is to show that with probability one, for each fixed s, 
(3.6) P{Hr, <Cnm + sm^/^) = F^ign s) + o(l) 

as n ^ oo. We will demonstrate these asymptotics by pointing out the necessary modifications 
to the argument in [GTWl]. 

All the —l/vj in our previous discussion are contained in the interval (— oo, ^], where ^ = 
1—1/6. (Recall that b is the maximum of the support of dF.) Let F„ be the empirical distribution 
function given by 

dFn = n"^ ^ Sp. 

and let ( • denote the integration with respect to Recall the Glivenko-Cantelli theorem, 
which says that, with probability one, F„ converges uniformly to F as n — oo. 

We first show that, under our standing assumptions, the quantities and ti„ of Lemma 3.1 
converge almost surely as n ^ oo to the corresponding quantities associated with the distribution 
function F. Recall that we set r = p/{l — p), p = r/{l + r). We remark that cq in the following 
lemma is the same as c in Theorem 1, and uq = {a — l)/a. The notation has changed to conform 
with (3.4) and (3.5), which are in turn chosen to connect with the saddle-point approach in 
[GTWl]. 



Lemma 3.2. The equation 

1 

a 



(l + ruo)V (no -1)2 
has a unique solution uq G (^, 0) and if cq is then defined by 

1 / ruo 

Co = a 



1 - Uq \ 1 + ^Uq 

we have Cq G (0, 1). 



Proof. The argument goes almost exactly as for Lemma 3.1. The assumption a < Oc is equiv- 
alent to a (r/(l + ru)^) > l/{u — 1)^ when u = ^, while a > a'^ yields the opposite inequality 
when ti = 0. □ 
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Note that one obtains u„ as uq, except that the expectation ( • ) is replaced by the expectation 



Lemma 3.3. Almost surely, Un — > uq and Cn ^ cq as n — oo. 



Proof. Integration by parts gives 



{l + rz) 



r'^ °' d r 



The derivative in the integrand is Tiniformly bounded for z in any compact subset of the comple- 
ment of (—00, ^]. Hence the expected value is continuous in F and differentiable for z ^ (—00, 
Moreover, 



d 



a 



dz\ \(l + rz)V 

is negative, hence nonzero, z = uq. The statement concerning Un therefore follows from the 
fact that Fn ^ F uniformly and the implicit function theorem. The assertion for c„ then follows 
by a similar integration by parts. □ 



Lemma 3.4. There exists a (deterministic) wedge W with vertex v > ^, bisected by the real 
axis to the left of v, such that, almost surely, the curves lie outside W for sufficiently large 
n. 



Proof. First we show that, if e is small enough, C„ is disjoint from the disc 

D{i, e) = {z: \z-i\<e}. 

Prom the facts that a'^{un) = a'^{un) = 0, a'^'{un) > 0, and a-'^{z) 7^ for 2; G (^, Un), it follows 
that an is strictly increasing in the interval (^, u„). Therefore, we can choose small enough 
£ > and (5 > so that (T„(^ + e) < cr„(u„) — 26 for all large enough n. In addition, if e is small 
enough, 

log \z-l\ + (c„ - 1) log |z| < log + £ - 1| + (c„ - 1) log \^ + s\ + S 

for all z G D{^, e). Now each |1 + rjz\, and so its logarithm, achieves its maximum on D(^, e) 
at the point z = ^ + e. By combining the last three observations, we see that everywhere on 
D{^, e) we have 

3fic7„(z) < cr„(^ + e) + S < (T„('U„) - S. 
Since 3ft a achieves its minimum on C~ at z = Un, the curve must be disjoint from -D(^, e). 
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Figure 5. Wedge W, angle ei, and disk D(^,e) as described in the proof of Lemma 3.4. 

For a small ei > (possibly much smaller than e), denote by W the wedge with vertex 
^ — e/2 bounded by the real axis to the left of — e/2 and the ray arg {z — ^ + e/2) = vr — ei. Our 
next step is to show that C~ is disjoint from W if ei is small enough. As Q an{z) is constant 
on the portion of C~ in the upper half-plane, 

n 

- arg (1 + rjz) + arg (z - 1) + (c„ - 1) arg z = c„7r, 
i=i 

where all arguments lie in [0, vr]. For z £ W, 

arg {z -1) + (cn - 1) arg z > Cn arg z > c„(7r - ei). 

Since b is in the support of dF, the strong law implies that, almost surely, at least a positive 
fraction rj of the —1/rj lie in the interval — e/2, ^] for large enough n. The contribution of 
these terms (and nonnegativity of the others) in the following sum provides a lower bound valid 
for z e W: 

n 

— > arg (1 + rjz) > arjiir — ei). 
n ^-^ 

Hence, for z gW', 

"^CTniz) > ar?(7r - Si) + Cn(vr - Si) = {arj + Cn) (vr - ei) > c^tt. 
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if £i is chosen to be small enough. Therefore, C~ is disjoint from the wedge W for small enough 
£i. By symmetry, C~ is disjoint from the reflection of W over the imaginary axis. We have 
shown that the curve is also disjoint from D{^, e) and the union of the disc and the two wedges 
contains a wedge of the form described in the statement of the lemma. 

This establishes the statement of the lemma concerning C~ . Since is to the "right" of 
C~ (it begins to the right and they cannot cross), the statement for C+ follows automatically. 
□ 



In the following lemma ctq denotes the function a associated with the distribution F, 
ao{z) = a (log (1 + rz)) + log {z - 1) + (cq - 1) log z, 

and 



go = \uo\ ^ |^f^o'(wo)| 



-1/3 



Lemma 3.5. Almost surely, zu'^{z) — > zu'q{z) uniformly outside the wedge W of Lemma 3.4- 



Proof. We have 

/TZ 
Y^p^d(F„(p) - Fip)) + c„ - Co 



= a {Fn{p) - F{p)) , ^ dr + Cn- cq. 

^0 



^6/(1-6) 

(l + rz)2 

The last term goes to by Lemma 3.3. The last factor in the integrand is uniformly bounded for 
r G (0, 6/ (1 — b)), z and z bounded. Thus z (t'^{z) — > z cfq{z) uniformly on bounded subsets 
of the complement of 1^. If z is sufficiently large and outside W then it is outside some wedge 
with vertex bisected by the negative real axis, and on the complement of any such wedge 

is uniformly bounded. Thus za'^{z) — > za'Q{z) uniformly throughout the complement of W . □ 



The preceding lemmas show that the curves are uniformly smooth, as we now argue. 
The function u'q{z) can have no other zero in the complement of (— oo, ^) than at z = uq. This 
follows from uniform convergence and the fact that the corresponding statement holds for the 
(Tniz). Thus the functions o''^{z) are uniformly bounded away from zero on compact subsets not 
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containing uq. As we move outward (i.e., away from u„) along C^, SJcr is constant and is 
increasing on C~ and decreasing on C+. It follows that if s measures arc length on the curves 
then, for z G C^, 

(o 7^ dz _ KM 

This shows that the are uniformly smooth on compact sets (to be more precise, the portions 
in the upper and lower half-planes are). Moreover, they are uniformly close on compact sets to 
the corresponding curves for the distribution function F. In particular, the length of C+ is 
0(1). 

To see what happens for large z on C~ observe that 

lim z u'^{z) = a + c„ > 

uniformly in n. This and (3.7) show that |z| is increasing as we move far enough out along C~. 
If r is an arc of C~ going from a to 6 then 



Wn{z)\ds = j^a'^{z)dz = anib) - cT„(a). 



Hence the length of T is at most |6 — a| times 

mincer K(z) I 

where [a, b] is the line segment joining a and b, as long as this segment does not meet (— oo, ^). 
It follows from the above, for example, that the norm of the function (1 + |2p)~^ on C~ is 
Oil). 

In [GTWl] we needed asymptotics with error bounds for all j,k < h and this required a more 
careful analysis of the integrals in (3.1) and (3.2) than we indicated; instead of the steepest 
descent curves passing through the same point they pass through different, but nearby, points. 
With what we now know we can show that these curves are uniformly smooth with uniformly 
regular behavior near infinity, and this is what is needed to see that in our case the asymptotics 
hold uniformly in n. 

Lemma 3.6. Almost surely, gn ^ 9o ^ '^''^d (3.6) holds. 

Proof. The first statement follows from Lemmas 3.3 and 3.5 and the fact that the c'^iz) have 
only two zeros outside (— oo, ^] counting multiplicity, and therefore a"'{uo) / 0. 
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To establish (3.6), one now has to go through the the steepest descent argument in Sec. 3.1.2 
of [GTWl], and make some obvious changes, justified by the results of this section. For the 
analogue of Lemma 3.1 there, for example, we would add the phrase "and all sufficiently large 

n" to the end of the statement. At the end of second sentence of the proof wc would add the 
phrase "since a'^{un) is uniformly bounded away from zero the length of is 0(1)." After 
the last sentence we would add "again since the length of is 0(1)." Analogous changes need 
to be made throughout the argument and we skip further details. □ 

Proof of Theorem 3. By Lemma 3.6, we can take G„ = c„m. □ 

Proof of Theorem 2. Note first that 

V 1 + Too y 

where Uq is as in Lemma 3.2 (and, as we remarked earlier, cq = c). 

The proof rests on the crucial property (3.6) and the fact that n^/^(F„ — F) converges in 
distribution to a Brownian bridge B with an appropriate covariance structure; in particular 
5 is a Gaussian random element in -D[0, 1] ([Bil], Th. 14.3). By the Skorohod representation 
theorem, we can couple and B on some probability space J7o so that 

(3.8) n^/'^{F^-F)^B 

in fact converges for every uj ^ ([Bil], Theorem 6.7). We now prove that, under this coupling, 
the solution of (3.4) satisfies 

(3.9) n„ = no + n-^/2C/ + o(n-^/2-)^ 
for every u and for some Gaussian random variable U. 

To establish (3.9), define 

1 

1 

~ {u-lf 

By Lemma 3.6 and its proof, there exists a (deterministic) neighborhood C C of no in which, 
with probability 1, n„ (resp. uo) is for large n the unique solution to 9n{u) = (rcsp. 6q{u) = 0). 
Therefore we can choose a fixed contour C in U such that and no are given by 

1 f , 1 f W . 

27^^ Jc On{u) 2m Jc 9o{u) 
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By (3.8), we have, uniformly for u E C, 



0„(u)=0o(«)+n-V2a(-^l— ) +o(n-V2). 

\(l + ra)2/i3 

Here ( • )g is the expectation with respect to dB, but by integration by parts (as in the proof of 
Lemma 3.3) we can make B appear in the integrand. Therefore 

dn{u) 9o[u) du 9o{u) 

If we multiply this identity by u/2Tri and integrate over C the left side becomes the first 
term on the right becomes uq while the second term becomes n~^^^ U where 

is a Gaussian random variable. This proves (3.9). 
Let 

1 / r 

^n{u) = z «\ T"^ 

1 — u \l + ru 

so that Cn = (pn{un)- We claim that 

(3.10) Cn = <fn{uo) + O{n-^). 

To see this, we use the fact that ip'^{un) = to write 

Cn = (fniUn) = V>n{uo) + {Un - Uof / t(p'^{t Un + {1 - t) Uq) dt. 

Jo 

Thus, (3.10) follows from (3.9) and the uniform boundedness of the (pniu) near u = uq. 
Now, by the central limit theorem. 



n 



iy 



_rjUo I ruo 
1 - 



converges in distribution to a Gaussian random variable X with mean and variance . There- 
fore, 

(3.11) Vn{cn - Co) aX. 
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Finally, (3.6) implies that, for any S > 0, 

(3.12) P{-Sm^/^ <H- Cr,m < Sm^/^) 1. 

(In fact, (3.6) implies that the above statement holds with probability 1 before the expectation 
with respect to pi, . . . ,p„ is taken, that is, if H is replaced by H^-) It follows from (3.11) and 
(3.12) that 

{H — Com) / y/m y/a X, 

which concludes the proof. □ 

Remark 1. We did not need the full force of (3.6) for the above proof to go through. Instead, 
a much weaker property (3.12) suffices. 

Remark 2. As mentioned in the Introduction, independence of p„ is not necessary for the results 
of this section to hold. Indeed, one only needs Glivenko-Cantelli theorem for convergence in 
probability of H/m to the time constant, hence ergodicity of pi,p2, • • • is enough. Furthermore, 
a strong enough mixing property of this sequence is sufficient for a normal fluctuation result. 
This follows from Billingsley's results in Section 22 of the first (1968) edition of [Bil]. 

Remark 3. We assumed that a = n/m is fixed, but the proof of Theorem 2 remains valid with 
n = am + o{y/m). 
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